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In this paper we develop a formalism for incorporating hadron loops in the quark model. We de- 
rive expressions for mass shifts, continuum components and mixing amplitudes of "quenched" quark 
model states due to hadron loops, as perturbation series in the valence-continuum coupling Hamilto- 
nian. We prove three general theorems regarding the effects of hadron loops, which show that given 
certain constraints on the external "bare" quark model states, the valence-continuum coupling, and 
the hadrons summed in the loops, the following results hold: (1) The loop mass shifts are identical 
for all states within a given N,L multiplet. (2) These states have the same total open-flavor decay 
widths. (3) Loop-induced valence configuration mixing vanishes provided that Li 7^ Lf or Si 7^ Sf. 
The charmonium system is used as a numerical case study, with the ^Po decay model providing the 
valence-continuum coupling. We evaluate the mass shifts and continuum mixing numerically for all 
IS, IP and 2S charmonium valence states due to loops of D, D*, Da and D* meson pairs. We find 
that the mass shifts are quite large, but are numerically similar for all the low-lying charmonium 
states, as suggested by the first theorem. Thus, loop mass shifts may have been "hidden" in the 
valence quark model by a change of parameters. The two-meson continuum components of the 
physical charmonium states are also found to be large, creating challenges for the interpretation of 
the constituent quark model. 

PACS numbers: 12.39.-x, 13.20.Gd, 13.25.Gv, 14.40.Gx 



I. INTRODUCTION 

The discovery of the narrow charm-strange mesons 
D*o(2317)+ [H and D,i(2460)+ Q has given special im- 
petus to the calculation of hadron loop effects, since the 
loops are often cited as a possible reason for the sur- 
prisingly low masses of these mesons. This possibility is 
supported by the prediction that the corresponding cs 
quark model states have especially strong couplings to 
the open-flavor decay channels DK and D*K [3, 4]. (For 
discussions of the importance of hadronic loops effects 
in this and other contexts, including other heavy-quark 



mesons J see for ex an^ le Ref s ■ [5|,^ | |8|, 1 1^ U^l 1 1 2| , 13, 

The subject of valence-continuum coupl ing s is also rel- 
evant to the X(3872) seen in J/ipTr~^Tr~ [33, [s^, which 
may be dominantly a weakly-bound l"*""*" DD* molecular 
state; the size of the cc valence component present in this 
system has important implications for the properties of 
this state. (See Ref. ^39j for a review of these and other 
recent developments in heavy-flavor hadrons.) 

Since the open-flavor decay couplings of hadrons to 
two-body final states A BC are large, one might antic- 
ipate that second-order decay loops, in which a hadron 
virtually decays to a two-body intermediate state and 
then reforms the original hadron (A — > BC — > A), are 
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also important effects. These second-order virtual pro- 
cesses give rise to mass shifts of the bare hadron states, 
and contribute continuum components to the physical 
hadron state vectors. A careful estimate of these mass 
shifts is of great interest, since they are usually not in- 
cluded in quark potential models and are only partially 
present in quenched lattice QCD, and may constitute im- 
portant "systematic" errors in the results. 

In our initial study we develop a formalism for treating 
these loops effects, using results from our earlier studies 
of open-flavor decay amplitudes. In particular we give 
results for the loop-induced mass shifts and continuum 
amplitudes of hadrons, as well as the off-diagonal "spec- 
troscopic mixing amplitudes" induced by hadron loops 
between different external discrete hadron basis states. 

As a numerical application we consider the charmo- 
nium system, and evaluate these mass shifts and contin- 
uum components for the lighter (narrow) charmonium 
states that lie below the open charm threshold. Charmo- 
nium is especially attractive as a test system for studying 
loop effects because the low-lying spectrum is clear ex- 
perimentally, with complete IS, IP and 2S multiplets, 
and all eight states in these multiplets are below the 
open-flavor decay threshold. (This implies that all mass 
shifts are negative, with no cancellations.) In addition 
the charmonium system is only moderately relativistic, 
and the spectrum is quite well described by quenched 
potential models and lattice gauge theory. Thus loop 
effects may be evaluated as (possibly) perturbative cor- 
rections to well understood cc potential model states, and 
the results may be contrasted with an unambiguous ex- 
perimental spectrum. 
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II. FORMALISM 



A. The Loop Model 



To incorporate hadron loop effects in tlie quark model, 
we model a physical hadron as a bare valence state |A) 
augmented by two-hadron continuum components, 



\^) = |A)+^^Bc|BC) 



(1) 



BC 



We assume that the Hamiltonian for this combined sys- 
tem consists of a valence Hamiltonian Hq (the quark 
model Hamiltonian) and an interaction Hi which couples 
the valence and continuum sectors, 



H = Ho + Hi 



(2) 



We will evaluate the continuum components of the 
hadron state and their physical effects as a perturbation 
series in the valence-continuum coupling Hi. Our start- 
ing point for this perturbations series is the set of single 
valence hadron Hq eigenstates; a specific valence state is 
written as |A(pa)), and is assigned an Hq eigenvalue of 
Ea = (M^ + p^y/"^. Since we normally work in the rest 
frame of the valence hadron, pA = 0, this energy eigen- 
value is just the rest mass Ma of the bare valence quark 
model hadron. 

The free two-hadron valence states which form our 
zeroth-order noninteracting continua are written as 
|B(pb)C(pc))- The valence Hamiltonian Hq is under- 
stood to operate only between the constituents of B and 
C separately; BC interactions, which are not treated here, 
would be incorporated in a separate two-hadron inter- 
action Hamiltonian. This BC continuum state has Hg 
eigenvalue Ebc = Eb -f Eq where Eb = (M| + piY^"^ 
and Ec = (M^ +p^y^'^. In the A rest frame we have 
Pb — —pc = P, and with p = \p\ the energies are 
Eb = {Ml + p2)i/2 and Ec = (M^ + p^f/^. 

The matrix elements of the valence-continuum cou- 
pling Hamiltonian are of the form 



(BC|Hi|A) -/i/,,5(pA-PB-Pc) 



(3) 



With explicit momentum labels these rest-frame one- 
and two-hadron valence states are written as |A(0)) and 
|B(p )C(— p )), and the coupling matrix element is a func- 
tion of a single momentum vector, hfi{p). 



B. Mass Shifts 

The mass shift of a valence hadron A due to its cou- 
pling to the BC continuum may be expressed in terms 
of the coupling matrix element hfi(p) of Eq.([3]) using 
second-order perturbation theory (for a general discus- 
sion see .40] ) . The usual discrete sum over inter- 
mediate states generalizes to a momentum-space integral 



over continuum states |B(p)C(— p)); the result for a sin- 
gle BC channel is 



AM 



(BC)_ 



V (En-Eo) 



<Fp 



(Ebc - Ma) 



Mb+Mc (Ebc 



pEbEc 



rfEBC pEbEc [ 2 
,c-Ma) Ebc 7'''^"^'' 



Ma 



(4) 



.Ma 



where V is the principal part integral. There is an im- 
plicit sum over any intermediate-state polarization labels 
in the squared Hamiltonian matrix element 

As a check of our central result Eq. (0]) , note that the 
imaginary part of the mass shift 



Im(AMr^) = -vr|^|dl7>,,| 



should be related to the total decay rate by 



Ebc=Ma 



r(A BC) 



-2Im(AM^^^^) 



(5) 



(6) 



The standard A ^ B -I- C decay rate formula given in 
Eq.(5) of Ref.[4l| is indeed consistent with this relation. 

If the initial hadron mass is below BC threshold (Ma < 
Mb -I- Mc) we do not encounter a singular energy de- 
nominator, and this mass shift is a real, negative definite 
integral over p, 



AM 



(BC) _ 



p^dp 



(Ebc - Ma) 



(7) 



If one considers mixing between the valence state A and 
several continuum BC channels, the total mass shift at 
this (leading) order in the valence-continuum coupling is 
the sum of the individual mass shifts due to each channel. 

Nonperturbative estimates of the mass shift can be 
made in the absence of final state interactions by sum- 
ming iterated bubble diagrams. The result is a full prop- 
agator of the form 



iG{s) 



(s - M2 - I](s)) 



(8) 



where S is the one particle irreducible self-energy of the 
meson in question. The propagator pole yields the me- 
son mass shift and width. Contact to our perturba- 
tive, nonrelativistic results can be made by identifying 
^r(s) = -Im(E(s)) and 2^5M{s) = Re(S(s)), and 
assuming that the width and mass shift are small rela- 
tive to the unperturbed meson mass. 

We have computed numerical pole positions with it- 
erated loops for the charmonium examples of the next 
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section using this formalism, and find rather smaU differ- 
ences in mass shifts relative to the single loop approxi- 
mation (typically AM(cc) changes by less than 10%). 

We have also examined the effect of mixing-induced 
coupling between states. Thus the denominator of the 
propagator becomes a matrix, {s — mf)dij — Tjij(s). Solv- 
ing this equation for the case of the J/ tp coupling to i/'' 
and tp" through DD, D^D^, DD*, D^D*, D*D*, and D*D* 
continua again yields corrections to the one-loop diago- 
nal results that are typically ca. 10%. We therefore sim- 
ply present perturbative (one-loop), single channel mass 
shifts in the discussion of charmonium. 



C. Continuum Components 

Although mass shifts due to loops may be "hidden" 
in fitted parameters in quenched approaches, such as mq 
or Vo in potential models and iriq or a(/3) in quenched 
LGT, it should nonetheless be possible to identify other, 
more characteristic effects of the two-meson continuum 
components. We will require the explicit continuum com- 
ponent wavefunctions to evaluate their effects on observ- 
ables. Here we give general results for these wavefunc- 
tions; an example will be considered in the discussion of 
charmonium. 

The valence-continuum coupling Hi induces a contin- 
uum component in an initially pure valence state |A). At 
leading order in hji this continuum component is given 
by 



5^ VbcIBC) = -(Ho - Ma)-^ Hi |A) 



(9) 



EC 



The momentum-space wavefunction of the continuum 
component in a specific channel BC is 



'/'B(p)C(-p) = 0BC(P) = - 



(Ebc(p) - Ma) 



(10) 



Using the conventions of Ref. [43| , the corresponding real- 
space wavefunction in the relative separation f = rs — ^'c 



(27r)3/2 



(11) 



For nonzero spin this spatial wavefunction is implic- 
itly summed over the meson orbital and spin magnetic 
quantum numbers, to give overall states with the J, Jz of 
meson A. 

The norm of this continuum component gives the prob- 
ability that the physical energy eigenstate is in the two- 
meson channel BC. This is 



(BC) _ |(V'n|gj|^0)P 



E 



(E„ - Eo)2 



^f^\ 



p^dp 



[Ebc - MaY 



(Ebc - Ma)2 ■ 

(12) 

dn^\hf,\'' . (13) 



D. Spectroscopic Mixing 

Mixing of discrete "valence" quark model basis states 
through hadron loops is an interesting effect which may 
have easily observable consequences. The amplitude a/i 
to find a discrete basis state |/) in the initially pure va- 
lence state \i) as a result of continuum mixing is given 
by second-order perturbation theory, 



1 



a/i 



(M; - M,) 



BC 



V,BC 



{p)hBcAp) 



(Ebc(p)-M,) 



(14) 



For an initial valence state within the continuum this is 
replaced by a principal part integral and the amplitude 
afi has an imaginary part, analogous to Eq[4l 

Note that this loop-induced mixing amplitude is some- 
what counterintuitive, in that it is nonsymmetric in gen- 
eral; 



(15) 



This disagrees with the simple picture of an orthogo- 
nal rotation between two basis states often used to de- 
scribe mixing in the quark model. (Examples include 
mixing between spin-singlet and spin-triplet axial vec- 
tor Ki and Di mesons, and between the |2'^Si) and 
pDi) charmonium basis states in the {■0'(3686), -0(3772)} 
system.) Since this is actually an infinite-dimensional 
Hilbert space rather than a two-dimensional one, it is 
of course not necessary that | a/i | = \aif\- Instead the 
valence state \i) that is closest to the continuum, and 
hence minimizes the valence-continuum energy denom- 
inator (Ebc — Mi) in Eq. p^ . will tend to experience 
the largest mixing. This will be illustrated in the next 
section. 



E. Three Loop Theorems 

In the Appendix we show that sums over sets of mesons 
within the loop under certain conditions gives very simple 
relations between the mass shifts, strong widths, and con- 
figuration mixing amplitudes due to hadron loops. Al- 
though these relations are not exactly satisfied in nature, 
they are sufficiently accurate to be relevant to realistic 
problems such as the charmonium examples we consider 
here. 

Provided that our conditions are satisfied, one may 
show that for the states {A} in a given Ni,Li multiplet: 

1. The mass shifts for all states {A} are equal. 

2. Their strong (open- flavor) total widths are equal. 

3. The configuration mixing amplitude afi between 
any two valence basis states i and / vanishes if 
Li ^ Lf or Si Sf. 

These conclusions hold to all orders if there are no final 
state interactions in the continuum channels. 
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To prove these loop results, we consider a sum over a 
finite set of intermediate (loop) mesons that runs over 
all mesons in a given N,L multiplet, taking on all al- 
lowed values of spin S and total angular momentum J. 
Examples_of_such loop sets_include '^SS"; { (DD), (DD*), 
(D*D), (D*D*) } and "SP"; { (DD^), (D*D5), (DDi), 
(D*Di), (DDi'), (D*Di'), (DD^), (D*DJ) }, where we 
have explicitly indicated antiparticles. The proof as- 
sumes that all members of the set of intermediate (loop) 
mesons which are summed over have the same mass, and 
for the first two "diagonal" conclusions (equality of mass 
shifts and total widths) we also require that the external 
mesons have a common bare mass and radial wavefunc- 
tion. In addition there are general conditions that the 
valence-continuum coupling must satisfy, which are dis- 
cussed in the Appendix. These constraints are satisfied 
by the '^Pq model that is used for illustration in this pa- 
per. 

The first conclusion suggests how these intrinsically 
large loop mass shifts can be hidden in the parameters 
of quenched models; since the largest effect of loops on 
the spectrum of states in the multiplets we consider is an 
overall downward mass shift in the multiplet c.o.g., this 
can be approximately parameterized through a change in 
the quark mass or through a constant in the potential. 

As as illustration of this theorem in a specific case, 
Table HIl shows the relative mass shifts of all the IP char- 
monium levels due to their '^Pq couplings to DD, DD* 
and D*D* meson loops, assuming that all the bare IP 
cc masses are identical, the D and D* meson masses are 
identical, and each flavor system has a common radial 
wavefunction. Although the individual continuum chan- 
nels DD, DD* and D*D* make different contributions 
to the mass shift of each meson, the summed mass shift 
from all three channels is identical for each of the four 
IP mesons. Thus, if the mesons are initially degenerate, 
they remain degenerate after these loop effects are in- 
cluded. This has been noted previously by Tornqvist fsS] ; 
the results presented here can be considered an elabora- 
tion of this original observation. Furthermore, Tornqvist 
and Zenczykowski have studied the analogous effect in 
baryons in Refs. [H, [s^. 

One may also see evidence for the no-loop-mixing re- 
sult for states with different L (conclusion 3 above) in 
our charmonium example. In Table IIIII we show the in- 
dividual DD, DD* and D*D* one-loop contributions to 
the mixing between pSi) |2'^Si) and pDi) charmonium 
valence basis states. Note that in some of the disfavored 
cases, such as pSi) — > pDi), there is an almost com- 
plete cancellation of the final pDi) amplitude, due to 
destructive interference between the DD, DD* and D*D* 
loops. This destructive interference between loops is still 
evident but less complete for mixing between the higher- 
lying states |2'^Si) and pDi), because they are quite close 
to DD threshold; this causes the energy denominators to 
vary widely between channels, so the mass constraints 
assumed in the theorem are strongly violated. 



III. NUMERICAL RESULTS: APPLICATION TO 
CHARMONIUM 

A. Mass Shifts 

To illustrate this formalism we will evaluate the effect 
of open-charm meson loops on the masses and compo- 
sitions of IS, IP and 2S charmonium states. We use 
the well established ^Pq model P, [3, El, [H, as 
the valence-continuum coupling Hamiltonian, and ne- 
glect two-meson interactions. Our general approach is 
very similar to an earlier study by Heikkila, Tornqvist 
and Ono [l^l , although we find somewhat larger loop ef- 
fects than reported by this reference. 

The ^Pq model treats strong decays as due to a bilin- 
ear quark-antiquark pair production interaction Hamilto- 
nian, Hi = which is normally evaluated 
using nonrelativistic quark model matrix elements. 

A diagrammatic technique for determining the valence- 
continuum coupling matrix element hfi between a me- 
son A and a two-meson state BC in the "^Pq model is 
given in Ref . [lil . We use this approach to determine the 
{hfi} A-BC valence- continuum matrix elements. Gaus- 
sian momentum-space quark model meson wavefunctions 
were used, with unequal light, strange and charm quark 
masses. For simplicity a common width parameter /3 
was assumed for all charmonium and open-charm me- 
son wavefunctions; tests of the overlaps of more realis- 
tic Coulomb plus linear plus smeared hyperfine wave- 
functions with Gaussians shows that this is a reasonable 
"zeroth-order" approximation. 



1. J/ip mass shifts 

As a first numerical example we consider the mass shift 
of an initial valence J/ip cc state mixing with the DD 
continuum. The hfi matrix element for the transition 
J/ip D(p ) D(— p) for a rest J/tp in polarization state 
m is given by 



23 1 
hf^ = — 



3r„ 



7 



33 l + r„ 7ri/4^i/2 



pe-pV3(i+.,.)Vi„(f7p) 



(16) 

where r„ — mn/rric is the light (n ~ u, d) to charm quark 
mass ratio, p — p/ (3, (3 is the simple harmonic oscillator 
(SHO) meson wavefunction width parameter (taken to 
be the same for all mesons in this work), and 7 is the 
dimensionless ^Pq pair production amplitude. On sub- 
stituting this hfi into the mass shift formula Eq.(I7]), and 
including a flavor factor of two for neutral and charged 
DD loops, we find 



AM 



(DD) 

i/i, 



36 V 1 + r„ ; 7rl/2 (2(p2 + ^2^)1/2 _ ^^^^) 



(17) 
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where /i = M//3 for each meson. 

Numerical evaluation of this integral using Mj/^ = 
3.097 GeV, Mq = 1.867 GeV, (3 = 0.5 GeV, 7 = 0.35 
(motivated by total widths; see Fig. 2 of Ref . |42| ) and 
r„ = mn/rric = 0.33/1.5 gives the result 



AM^/°^ = -23.1 MeV. 



(18) 



Using the experimental J mass as the input bare mass 
in this manner is of course only appropriate as an esti- 
mate of the size of these effects. Since this is in effect a 
renormalization problem, the sum of the (unobservable) 
bare mass and mass shift should be identified with the 
experimental J /ijj mass. 

Although this DD contribution is a relatively small ef- 
fect, incorporation of higher (1S)(1S) channels shows that 
the summed loop mass shifts are quite large. The formu- 
las for the DD, DD* and D*D* loop integrals in the J/^ 
system are identical, but the relative spin-flavor factors 
of 1:4:7 give a combined mass shift that is an order of 
magnitude larger than for the DD channel alone. (These 
1:4:7 spin-flavor factors were reported earlier by Heikkila 
et al. QT^ for loop contributions to mass shifts, and by 
De Rujula et al. and Close [1^ for charm production 
cross sections.) On including all six D, D*, D^ and D* 
pair channels (with Md = 1.867 GeV, Md- = 2.008 GeV, 
Md, = 1.968 GeV, Mdj = 2.112 GeV and = 0.55/1.5), 
we find 



E 

11=1 



-457.5 MeV, 



(19) 



This very large mass shift appears to invalidate the 
quenched quark model. In the next section we will see 
that this scale of mass shift is actually common to all 
the low-lying charmonium states, and can therefore be 
approximately subsumed in a change of parameters (such 
as the charm quark mass rric or an overall constant Vq in 
the cc potential) . 



2. Mass shifts of other charmonium states 

One can understand how such large mass shifts may 
have been accommodated in pure cc quark models by 
evaluating the mass shifts of the remaining low-lying 
charmonium states below DD threshold. We again set 
the bare masses equal to the experimental values to gen- 
erate this estimate; the values used are M^/ = 3.686 GeV, 

M 



3.637 GeV, M^^ = 3.556 GeV, M^, = 3.511 GeV, 
- 3.415 GeV, Mh, = 3.526 GeV, M^^ = 2.979 GeV, 
and the other model parameters are as before. 

The resulting mass shifts are given in Tabled and ev- 
idently are all quite large. Note however that they are 
rather similar, so there is a much smaller scatter about 
the mean shift; the mean and variance are respectively 
-471 MeV and 49 MeV. The scatter of mass shifts within 
a multiplet is even smaller; the variance within the IP 



multiplet for example is just 24 MeV. (The similarity of 
mass shifts within a multiplet was discussed in the pre- 
vious section, and is a consequence of the general nature 
of the valence-continuum coupling model.) 

The large overall shift could be parameterized in a pure 
cc "quenched" potential model through a shift in rric or 
through the addition of a large negative constant Vq to 
the cc potential. One expects that the goodness of fit to 
the cc spectrum is rather insensitive to these modifica- 
tions. 

The J/ip — rjc and ?/;' — 77^ loop-induced mass splitting 
has been discussed previously by Eichten et al. 1} . These 
authors sum over the same set of intermediate states em- 
ployed here, but use the Cornell decay model for the 
strong decay interaction. They find a small loop-induced 
J/ip — r]c mass splitting of -3.7 MeV and a ■0' — f?c splitting 
of -20.9 MeV, bringing their model into good agreement 
with the experimental ijj' — rjc mass difference. Table U 
shows that we find a numerically similar ?/;' — rj'^ splitting 
of -24 MeV; however, the ground state mass difference 
due to coupling to the continuum is -34 MeV, indicating 
that loop effects induce a larger -0' — rj'^ mass difference 
of approximately -t-10 MeV. This is consistent with the 
bare model employed here, which finds a substantially 
smaller bare t/i' — 77^ mass differenceji^] than that of Ref. 



B. Continuum Components 

Although the large negative mass shifts may be "hid- 
den" by the choice of rric or Vq in potential models and rric 
or a{(3) in quenched LGT, it should nonetheless be possi- 
ble to identify other observable effects of the two-meson 
continuum components, since according to Table 1 their 
occupation probabilities are comparable to the valence 
cc components. To illustrate this we will evaluate some 
of these continuum component wavefunctions explicitly, 
and consider their effect on some experimentally observed 
properties of charmonium states. 

Recall from Eq llOl that the continuum component 
wavefunction in momentum space, 4>bc{'p)^ is given by 



(P ) = - 



hfi 



(Ebc(p)-Ma) 



(20) 



Again specializing to the DD component of the J/ij] as 
our example, this momentum space wavefunction is 



0dd(p) = '/'DD(p)Yi,„(ilp) 



(21) 



where 
8 

^ 27 



l + 3r^ 
l + r„ 



pe 



7rl/4^3/2 (2(M2,+p2)l/2_M^/^) 



(22) 
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TABLE I: Mass shifts (in MeV) and cc probabilities for low-lying charmonium states due to couplings to two-meson continua. 
This one-loop estimate sets the unperturbed bare masses to the experimental values, and assumes "^Po model and SHO wave- 
function parameters 7 — 0.35 and l3 — 0.5 GeV and quark mass ratios r„ — nin/nic = 0.33/1.5 and = ms/rric — 0.55/1.5. 



Bare cc State 




Mass Shifts by Channel, AMi 


(MeV) 






Multiplet State 


DD 


DD* 


D*D* 


D,D, 


d,d: 


d*d: 


Total 


^ cc 


IS J/i/'(l^Si) 


-23 


-83 


-132 


-21 


-76 


-123 


-457 


0.69 


r?c(l'So) 





-114 


-105 





-106 


-98 


-423 


0.73 


2S V'(2''Si) 


-27 


-84 


-126 


-19 


-70 


-113 


-440 


0.51 







-118 


-103 





-102 


-94 


-416 


0.61 


IP X2(1^P2) 


-40 


-105 


-144 


-33 


-88 


-111 


-521 


0.49 


Xi(l'Pi) 





-127 


-148 





-90 


-130 


-496 


0.52 


Xo(l'Po) 


-57 





-196 


-34 





-172 


-459 


0.58 


/lc(l'Pl) 





-149 


-130 





-118 


-107 


-504 


0.52 



TABLE IL Relative one-loop mass shifts of IP charmonium states in the equal mass limit. 



Bare cc State 



Relative Mass Shifts, AMi(LBc)/AMtot(LBc) 





Lbc = 






Lbc = 2 






DD DD* 


D*D* 


DD 


DD* 


D*D* 


1^P2 





1 


3/20 


9/20 


2/5 


l^Pl 


1 








1/4 


3/4 


I'Po 


3/4 


1/4 








1 


l^Pl 


1/2 


1/2 





1/2 


1/2 



TABLE III: Valence configuration mixing amplitudes afi due to loops (DD; DD*; D*D* ^Pi; D*D* ^Pi) in the 
{|^Si),|2^Si),|^Di)} system. The total Ufi is the sum of the individual loop contributions, as indicated. The labels |I) 
etc. refer to the physical (unnormalized) states one finds due to loop-induced mixing between |cc) valence states. (These are 
perturbative, one-loop results, with parameters as in Table U]) Note the approximate cancellations in AL 7^ mixing, and the 
non-symmetric mixing amplitudes. 





|I> 


III) 


|III> 






[1] 




-.013- 


.011 -i 


- .000 -i 


- .006 


= -.018 


-.089 - 


.017i + .086 - 


.010 - 


.020 = 


-.033 - 


.017j 




-.003 - 


.014 - .001 - 


- .026 = -.045 






[1] 






-.572 - 


.138i + .573 - 


.072 - 


.143 = 


-.214 - 


.138i 




-f.015 - 


.026 + .004 - 


h .008 = +.001 


+.340 - 


.469 4 


- .063 4 


- .126 


= +.060 






[1] 
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Note that this component formaUy diverges as Mj/^ 
2Md, due to a vanishing energy denominator; this shows 
that as expected the largest continuum components arise 
in valence states that are closest to the continuum. The 
spatial wavefunction corresponding to (/)dd(p) is shown 
in Fig. [11 



the interpretation of the constituent quark model and 
quenched QCD, which both neglect meson loops. The 
main issue is whether such large loop effects can be ab- 
sorbed into parameter redefinitions when computing ob- 
servables. 



D. Spectroscopic Mixing 




0.2 0.4 0.6 0.8 1 1.2 1.4 1.6 1.8 2 
r(fm) 

FIG. 1: DD Continuum Component Wavefunction of the J/^. 



As noted previously, discrete charmonium levels be- 
low the continuum mix at second order in the valence- 
continuum Hamiltonian Hi through hadron loops, pro- 
vided that both the initial and final valence states \i) 
and I/) have nonzero matrix elements to at least one 
continuum intermediate state |BC). 

Here we shall illustrate this effect by calculating the 
amount of mixing between low-lying 1 states. First 
we consider the pSi) and |2'^Si) cc basis states, which 
at leading order are identified with the J/4> and '(/'(SBSG) 
respectively. We will give explicit formulas for mixing 
through DD intermediate states, and simply quote nu- 
merical results for mixing through higher two-meson con- 
tinua. 

The hfi matrix elements required to evaluate these 
mixing amplitudes are 

hBC.CSi ^ DD) - 



C. J/tp continuum probabilities 

The probability of finding the physical J/ip in the DD 
continuum (from Eq. (jl3p ) is 

p(DD) _ 



2^ 

36 



ri/2 



(2(p2+;,2^)l/2„^^/^)^ 



(23) 



0.021 



Although this appears to be a reassuring small correction 
to the valence quark model description of the J/ijj as a 
pure cc state, when we calculate the probability that the 
physical state is in any of the i = 1 ... 6 meson contin- 
uum states DD, D*D, D*D*, D^D^, D^D/, D/D/, we 
again find that the summed contribution is quite large. 
Expressed as the probability that the physical J/4' is in 
the valence cc state, we find 



,(cc) 

J/lp 



(i) 
J/V 



0.685 



(24) 



i=l 



Just as was the case for the mass shifts, we find that 
the continuum components of charmonium states are 
very large. This represents an interesting challenge in 



23 1 + 3r„ 7 



33 l + Tn 7ri/4;3i/ 



/i/,Bc(DD ^ 23Si) 



r_ p^-p /3(i+.„) Y,^(^^) ^ (25) 



25/2 
3572 



^ 2 1 8 (1 + 3r„) 2' 

9(l + r„) 27(H-r„)3^ 



and 

;i/3c(DD-.3Di) = 

211/251/2 



39/2 



r 2 (l + 3r„) , 



(1 + r) 15(l + r„)3 



7^1/4^1/2 P'^ ^IraK^h) ■ l^U 

Substitution of these expressions in Eg. lfH)) and eval- 
uation of the overlap integral gives the l'^Si-23Si mixing 
amplitudes a/^. Our numerical results, using the same 
parameters and masses as previously, are given in Ta- 
ble ITTll 
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IV. SUMMARY AND CONCLUSIONS 

In this paper we presented a formahsm for "unquench- 
ing the quark model" through the incorporation of the 
effects of hadron loops on valence quark model states. 
We gave expressions for the mass shift, continuum com- 
ponents of the hadron state vector, and mixing ampli- 
tudes between discrete valence states that follow from 
hadron loop effects for a given valence-continuum cou- 
pling Hamiltonian. 

As a numerical example we applied this formalism 
to the experimentally well-established light charmonium 
system, using the ^Pq decay model for the valence- 
continuum coupling. We evaluated the mass shifts and 
composition of the physical charmonium states for all 
IS, IP and 2S states using perturbation theory in the 
valence-continuum coupling; these mass shifts and two- 
meson components were found to be quite large. Since 
the mass shifts of the different charmonium levels are 
numerically rather similar, we speculate that they have 
been hidden in the choice of rric or a constant potential 
shift Vb in cc valence potential models. It remains to be 
seen whether the two-meson continuum components can 
be "parametrized away" - it is possible that they lead to 
important mixing effects between discrete charmonium 
basis states that may be experimentally observable. 

The mixing effects we find using the ^Po decay model 
as the valence-continuum coupling prove to be quite 
large for higher-mass intermediate continuum states. Al- 
though it is possible that these effects can be largely 
renormalized away, an accurate description of loop ef- 
fects will probably require the development of a more re- 
alistic valence-continuum coupling Hamiltonian than the 
^Pq model. 



APPENDIX: LOOP THEOREMS 

Numerical experiments suggest that although individ- 
ual loop contributions to physical observables are large, 
in practice there are often important cancellations or con- 
straints when loop sums over sets of mesons are carried 
out. This is evident for example in the mass shifts in 
Table HI the individual loop mass shift for a given state 
varies widely between states, but the total mass shifts 
when summed over loops are rather similar. One can see 
that these relations are exact in certain limits. As an ex- 
ample. Table HTl shows the relative mass shifts of the four 
P-wave charmonium states in the limit in which they 
have identical initial masses, and the D and D* within 
the loops also have identical masses; although the indi- 
vidual channel mass shifts differ, we find the same total 
mass shift for each P-wave state on summing over the 
channels DD, DD* and D*D*. 

A similar result is evident in the loop-induced config- 
uration mixing discussed in the text; the configuration 
mixing amplitude a fi between initial i and final / meson 
basis states in the usual N, J, L, S basis is found to be zero 



if Li 7^ Lf, provided that the mesons in the loops have 
identical masses and we again sum over a complete set 
of loop meson spin states Sb and Sc. As an example, in 
this limit this gives a zero mixing amplitude due to loops 
between any charmonium ^Si and "^Di basis states. 

In this appendix we give a proof of this mass shift 
identity and the zero-mixing result for loop sums; these 
results hold whenever one sums over loops containing 
a complete set of spin (S) meson states (in a given 
N,J,L,S multiplet). The proof applies to the "^Po cou- 
pling model in particular, but also holds for a more gen- 
eral class of valence-continuum couplings, specifically to 
spin-one, factorised, spectator decay models, as discussed 
by Burns, Close, and Thomas fsO*!. In this type of model 
the valence-continuum coupling proceeds through spin- 
one qq pair production, the initial quarks do not couple to 
the decay vertex, and the spatial dependence of the decay 
vertex multiplies the created qq spin operator: O — aTp, 
where ip represents the spatial portion of the decay ver- 
tex. The proof therefore also applies to the Cornell decay 
model [? ] and a decay model based on the nonrelativis- 
tic reduction of the interaction J 'il)'il}{x)V{x — y)'il)ip{y), 
but does not apply to pair production from one gluon 
exchange (discussed in Ref.|4l|). 

Given a valence-continuum coupling of this general 
form, which includes the "^Pq used in this paper for 
numerical examples, one may show that the general 
(BC|IIi|A) matrix element is of the form 

{Ja [LJbc] ; JBC [jB jc] ; JB [sb^b] jc [sc^c] ki/> I Ja [sa-^a]) = 

(-)''iL/SBC^BciBjcjBC5ASBSC5BC ' 

(%[L£bc];^bcMc]||^II^a> • 

SB ^B Jb ] ( 1/2 1/2 SB ] 

sc JG > < 1/2 1/2 sc \ ■ 

SBC ^BC jBC } { SA 1 SBC J 
/ SBC ^BC jBC \ / SBC SA 1 \ ,^ 

\ L JA Lf j\ £a Lf M / ^'^■^> 

where x — ^/2x + 1 and r] = L + sbc + ^bc + Lf -t- sb . 

If the expressions for the mass splitting (EqlT]) or spec- 
troscopic mixing are summed over intermediate states 
BC with identical masses, the resulting common energy 
denominators may be taken outside the sum over chan- 
nels, and one is left with the expressions 

J (2ny [m, - Ebc(p) + «e) ^ 

(A.2) 

and 

_ 1 r (fp 1 

" (m,-m/) J (2^)3 (TO,-EBc(p)+ie) ' 

Y,{jf[sfef]\<TiP\BC){BC\aiPUsie,]) (A.3) 

BC 



9 



where Ebc is the common energy of ah states in the same 
multiplet as BC. 

The sum over intermediate states simphfies when one 
considers a subsum over spin muhiplets: 



(A.4) 



BC 



sbscJbJc 



the angular momenta Ib, L can remain fixed. 

On substituting Eg lA. II into Eqs lA.2l and IA.31 and us- 
ing the orthogonahty relation for 9j and 6j symbols 



Jl J2 Jl2 

h ji hi 

jl3 j2i J 



E •?13-^24 

il3 j"24 




.?12.?34 



and 



E^- 



J12 





32 


jl2 \ 




32 


il2 \ 


\J3 


J4 


J J 


\J3 


J4 


J' i 



5{J, J') 
J ' 



we obtain the following sum: 



(A.5) 



(A.6) 



(A.7) 



Since this expression is independent of the initial and 
final meson spin, we conclude that all mesons in a given 
(assumed degenerate) spin multiplet receive the same 
width and mass shift from the sum over all intermediate 
(loop) mesons in a given spin multiplet. Furthermore, 
the spectroscopic mixing between mesons of different or- 
bital angular momentum is zero when sums over spin 
multiplet intermediate states are carried out. (The ex- 
ternal meson masses need not be identical to prove this 
result.) Finally, since these matrix elements drive non- 
perturbative mixing (see the discussion following Eq|5]), 
these conclusions also apply to nonperturbative mixing, 
in the absence of final state interactions. 

Spectroscopic mixing between mesons with differing 
radial quantum numbers (but identical otherwise) is not 
zero in general. The size of this mixing is governed by the 
spatial dependence of the strong decay vertex. Spectro- 
scopic mixing has been studied previously by Geiger and 
Isgur.29j, who considered the closure approximation, in 



which all loop mesons are assumed to be degenerate, not 
simply those in a spin multiplet. Geiger and Isgur used 
this approximation to explain the observed weakness of 
loop-driven OZI violation effects. We remark that the 
closure approximation implies that spectroscopic mixing 
between states with different radial quantum numbers is 
zero: under this approximation Eq lA. 31 simplifies to 



1 



(nii - ruf) {m, - E) 

(O|02|O) 
(to; - mf){mi - E) 



{nfjf[sf£f]\0'\n,J,[sd^) 



{nfjf[sfef]\ndi[si£i]) 



(A., 



where the last form follows from the spectator nature of 
the decay model. Thus, if we impose the equality of all 
loop meson masses, all spectroscopic mixing is zero in the 
■^Pq model, and in a wide range of related decay models. 

Finally, the previous discussion remains largely un- 
changed when considering mixing between initially de- 
generate states. In this case one must diagonalize the 
matrix of second order matrix elements in the degen- 
erate subspacefsH. SHij — {rrii — mj)aji. Under the 
conditions of the theorem, off-diagonal matrix element 
in SHij are zero when the meson spins or angular mo- 
menta differ. Furthermore the diagonal matrix elements 
are identical. Thus conclusions concerning mass shifts, 
widths, and small or zero spectroscopic mixing remain 
unchanged. 
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